The ground state of a class of noncritical ID quantum spin systems can be 

approximated efficiently 
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We study families H n of ID quantum spin systems, where n is the number of spins, which have a 
spectral gap AE between the ground-state and first-excited state energy that scales, asymptotically, 
as a constant in n. We show that if the ground state |fi m ) of the hamiltonian H m on m spins, where 
m is an O(l) constant, is locally the same as the ground state |fin), for arbitrarily large n, then 
an arbitrarily good approximation to the ground state of H n can be stored efficiently for all n. 
We formulate a conjecture that, if true, would imply our result applies to all noncritical ID spin 
systems. We also include an appendix on quasi-adiabatic evolutions. 

PACS numbers: 75.10.Pq, 03.67.-a, 75.40.Mg 



I. INTRODUCTION 

The physics of low-dimensional lattices of quantum 
spins is especially rich and varied. As a consequence, a 
great deal of effort has gone into understanding the stat- 
ics and dynamics of these systems. However, despite this 
effort, many fundamental questions about quantum lat- 
tice models remain unanswered. Perhaps one of the most 
pressing of these questions is: can we faithfully approxi- 
mate, at least via some efficient numerical procedure, the 
ground-state properties of quantum lattice systems? If 
this were not true for physically realistic models then we 
would have to give up any hope of extracting theoretical 
predictions from these models. 

It is a folk theorem that, at least in one dimension, an 
approximation to the ground-state properties of a non- 
critical [2c| chain of quantum spins may be obtained ef- 
ficiently on a classical computer. Significant progress to- 
wards proving this theorem has been obtained recently in 
a ground-breaking paper by Hastings d| . Hastings found 
a procedure whereby an approximation to the ground 
state of a noncritical quantum spin system could be ob- 
tained and stored using subexponential resources that 
scale as n clog( ™) , where c is some constant which depends 
on the spectral gap AE and the local spin dimension. 
The computational complexity of this method is not far 
off the expected result, i.e., n r , where r is some constant. 

There is one procedure which appears to provide ar- 
bitrarily good approximations to the ground-state prop- 
erties of ID noncritical quantum spin systems, namely, 
the density matrix renormalisation group (DMRG). (See 
and references therein for a description of the DMRG 
and its relatives.) The DMRG is the premier tool used 
in numerical explorations of the physics of ID quantum 
systems, and has been used with unparalleled success in 
simulating both their statics, and more recently, dynam- 
ics. Many exciting extensions of the DMRG have been 
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developed, including, a powerful variant for 2D systems 
I. 

Unfortunately the DMRG is not known to be correct. 
That is, it is unclear if the DMRG always faithfully re- 
turns an approximation to the ground state and not some 
other low-lying excited state. Additionally, the complex- 
ity of the DMRG is currently unknown. It is entirely 
possible, in principle, that the DMRG requires exponen- 
tial resources to obtain a faithful approximation to the 
ground-state of noncritical spin system. However, exten- 
sive numerical experimentation strongly indicates that 
the DMRG requires only linear resources in n to sim- 
ulate noncritical systems. Nevertheless, hard instances 
for variants of the DMRG do exist Q , which means that 
we must be cautious when applying the DMRG in certain 
situations. This strongly suggests that while the average- 
case complexity of the DMRG may be polynomial, the 
worst-case complexity is probably exponential. 

The DMRG can be thought of as a variational min- 
imisation of the energy over the class of finitely corre- 
lated states or matrix product states (for an introduction 
to finitely correlated states and a detailed description of 
their properties see [f|). Finitely correlated states are 
particularly well-suited to this task because: (a) it is 
expected that they approximate realistic ground states 
well; and (b) there is an efficient computational proce- 
dure to extract local properties, like correlators, from a 
state stored as a finitely correlated state. The validity 
of the FCS ground-state ansatz is conditioned, at least, 
on the truth of (a), thus it is very desirable to show that 
good approximations to the ground states of some phys- 
ically interesting class of spin systems could be stored 
efficiently as a FCS. 

There are many ways to obtain an approximation to 
the ground state of a quantum system. For example, in 
the case of the DMRG, there are variants @, 0, [1] which 
obtain ground-state approximations via imaginary time 
evolution. However, in this paper, we'd like to emphasise 
another method to obtain ground-state approximations, 
namely, via adiabatic continuation. The idea with adi- 
abatic continuation is to start with a hamiltonian H(0) 
whose ground-state is known exactly, and then to adia- 
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batically vary along a path of hamiltonians H(s) until 
the desired hamiltonian H(l) is reached. Under the adi- 
abatic dynamics, the state at the end of the evolution 
will be the ground state of the final hamiltonian H(l). 

The spectral gap AE(s) between the ground- and first- 
excited state energy of H(s) provides the fundamental 
obstruction to approximating adiabatic dynamics: the 
smaller AE(s) is, the harder it is to approximate the 
dynamics. At this point we'd like to point out an ob- 
vious (but crucial) fact: if arbitrary paths H(s) are al- 
lowed then AE(s) can be made as large as desired up to 
min{Ai5(0), AE(1)}. However, in the context of spin 
systems, we don't allow arbitrary paths because they 
would presumably lead to an unphysical situation where 
H(s) contains interactions between many spins. Rather, 
we demand that H (s) has only local interactions through- 
out the path s £ [0,1]. This additional constraint mo- 
tivates us to define the notion of adiabatic connectivity: 
two quantum spin systems H and K are said to be adi- 
abatically connected if there exists a path of local hamil- 
tonians L(s) such that L(0) = H and L(l) = K, and the 
spectral gap for L(s) satisfies AE(s) > for all s G [0, 1]. 

In the case that a hamiltonian H is adiabatically con- 
nected to another hamiltonian K via a path L(s) with 
AE(s) > const, it turns out that ground-state proper- 
ties of H can be efficiently and certifiably obtained from 
those of K 0| . Because of certain counterexample sys- 
tems for DMRG methods [29[ it appears that adiabatic 
continuation is the only method whereby certifiable ap- 
proximations to the ground-state can be obtained effi- 
ciently. Thus, the problem of understanding the ground- 
state properties of a quantum spin system H can be re- 
duced to finding a hamiltonian for a well-understood spin 
system K which is adiabatically connected to H . 

In this paper we consider the problem of proving that 
the isolated eigenstates of a certain class of noncritical 
quantum spin systems can be efficiently represented as 
finitely correlated states with polynomial computational 
storage resources (in n and 1/AE). (The reason we 
say "isolated eigenstates" here is because our subsequent 
derivations make no use of the fact that the eigenstate 
in question is the ground state. For example, the argu- 
ment applies equally to the highest-energy eigenstate.) 
The noncritical systems we consider are local hamilto- 
nians which satisfy a crucial additional requirement: we 
assume that the ground state |O m ) of the system on m 
spins is locally close to the ground state |fi n ) for n spins, 
where n is arbitrarily large. This is a fundamental physi- 
cal assumption which, philosophically, underlies the suc- 
cess of the DMRG and relatives. In the case that this 
requirement is satisfied we show H n is adiabatically con- 
nected to a hamiltonian K whose ground state is exactly 
and efficiently known. As a consequence, if we know the 
ground-state energy f2„ for all n, we show that the ground 
state of H n may be efficiently approximated by a finitely 
correlated state. 



II. FORMULATION 

We will, for the sake of clarity, describe our results 
mainly for a finite chain C of n distinguishable spin-1/2 
particles. The family H of local hamiltonians we focus 
on (which implicitly depends on n) is defined by H = 
Y^j=o hj i where- hj acts nontrivially only on spins j and 
j + 1. We set the energy scale by assuming that \\hj\\ 
scales as a constant with n for all j = 0, 1, . . . , n — 1, 
where || • || denotes the operator norm. We can easily 
accommodate next-nearest neighbour interactions etc. by 
blocking sites and thinking of the blocks as new (larger) 
spins. However this can only be done a constant number 
of times: the quality of our approximation will decrease 
exponentially with the number of such blockings. We do 
not assume translational invariance. 

We make three major assumptions about our system. 
The first is that the spectrum of H = Ylj=o Ej\Ej)(Ej\ 
has a spectral gap AE = E\ — Eg between the ground- 
state energy and the first-excited state energy which is 
always strictly positive and scales as a constant with n. 
(This is the noncriticality assumption.) The second as- 
sumption we make is that the ground-state energy is set 
to zero. This is the principle reason why our analysis 
does not allow us to, in principle, efficiently calculate ap- 
proximations to ground-state properties of H because to 
perform this operation we need to know Eq — it is po- 
tentially a computationally difficult task to approximate 
the ground-state energy eigenvalue [l(| HH, 03 • Our final 
assumption is that the ground state |fi m ) for m spins, 
with m an 0(1) constant, is locally similar to \Q n ), with 
n arbitrarily large. This means that there exist unitary 
operators U and V which act nontrivially only on a con- 
tiguous block Ai (respectively, A2) of a constant number 
I <C m of spins located at the left (respectively, right) 
end of the block of m spins such that fidelity 

X = (n m \tf ®i®vip£ ) u®i®v\n m ) (1) 

is an O(l) constant independent of n, where p"m = 
tiffi(\Q, n ) (ft n \) is the reduced density operator for \fl n ) 
on m contiguous spins. In words: we assume that the 
ground state |f£ m ) of H m has some overlap with |f2„) 
when we are allowed to apply some correction operations 
to the ends of the chain of m spins. The physical idea un- 
derlying this assumption is that for noncritical spin sys- 
tems the ground state of m spins ought to be the same as 
that for n > m spins apart from boundary effects which 
should persist only a distance / = c/AE, with c some con- 
stant, into the bulk of the ground state of both systems. 
While this is physically reasonable we've been unable to 
show that it's true for all noncritical spin systems. Hence 
we make this an assumption. 

We are going to make a further simplifying assump- 
tion about the systems we are considering, namely that 
the overlap x' = {Q m \p"m \Sl m ) is an O(l) constant. 
This is obviously a far stronger assumption than that 
of the Eq. Jl}. However, it turns out that this as- 
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sumption entails no loss of generality in our subsequent 
derivations. The way to see this is to first notice that 
the state U ® I ® V\fl m ), where U and V are chosen 
as in Eq. (JTJ), is the unique gapped ground state of 
H' m = U <g> I <g) VH m W <g> I (g) 0. As long as f/ and 
V act on only a small number of spins (in comparison 
to to) near the boundary then H' m will also be a local 
hamiltonian. Our subsequent analysis only requires that 
our start hamiltonian H' m contains local interactions. 

Before we end this section we introduce some notation 
for approximations. If we have two quantities A and B 
then we use the notation A < B to denote the estimate 
A < CB for some constant C independent of n. Because 
we'll be interested in the consequences of allowing the 
minimum gap AE and the overlap x to depend on n 
we'll explicitly retain any dependence on AE and x in 
our calculations. 
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Adiabatic dynamics (U): 
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\Cl 2m ) = U\Cl m )\Q m ) 



FIG. 1: Illustration of the procedure to glue two copies of the 
III. ADIABATIC CONNECTIONS BETWEEN ground state \Q m ) for H m together via adiabatic continuation 

H m AND H„ to obtain an approximation to the ground state of H2m- 



In this section we construct a hamiltonian K, whose 
ground state is known exactly and which is adiabatically 
connected to H n . We construct K iteratively: we show 
that the hamiltonian K m , which we define to consist of 
two copies A and B of H m , is adiabatically connected 
to Hi m . Thus we "glue" the ground states of A and B 
together via adiabatic continuation. We then show that 
the adiabatic continuation from K m to H2 m can be ap- 
proximated by a unitary operator which acts on only a 
constant number of sites across the boundary between 
A and B. We then iterate this gluing procedure to ob- 
tain the ground state of H n . See Fig. [1] for a schematic 
illustration of our procedure. 

We describe our method in a slightly more general con- 
text: we show how to glue together two systems A and 
B of to and 77 — m spins, respectively. The way we show 
that a good approximation \Q) to the ground state \Q) 
can be stored efficiently is to consider the (adiabatic) dy- 
namics of an auxiliary system K. which is constructed in 
the following way. First fix n, the total number of spins. 
Next partition the chain into two contiguous blocks A 
and B. Now consider the following hamiltonian 

K = H — Hj (2) 

where Hi is an interaction term, to be defined, which 
spans the boundary between the two blocks. (Hj will not 
be the same as hi, the interaction term in H which spans 
the boundary.) If A consists of the first to spins, m < n, 
then B consists of the last 77—777 spins and we define Hi = 
hi + 51, where 8 is some constant. Thus we can write 
K = Ha + Hb — SI, where Ha — Y^j=a hj and Hb = 
Y^jZm ■ ^ e note that K has a unique ground state 
and spectral gap AE because of the assumed spectral 
structure of the family H n , i.e., both Ha and Hb belong 
to the family H n (Ha ~ H m and Hb ~ i? n _ m ). We set 



the constant 6 so that the ground-state energy of K is 
zero. (Recall that we've already set the zero of energy by 
requiring the ground-state energy of H is 0.) 

Now we construct a new system C whose hilbert space 
is a direct sum of two copies of the old hilbert space. The 
hamiltonian L for the new system is a direct sum of H 
and K : 

(3) 

The hilbert space for our new system is thus given by 
He = ®?=o W e think of this hilbert space as that 
of the original chain C of 77 spins with an extra spin, which 
we call C , that lives between spins to — 1 and to. Thus 
we can write L as L — Ic> <8> K + ^c+^c" j ^ jj j wnere 

We next observe, by the assumed properties of H 
and K, that the spectrum sp(i) of L has the follow- 
ing structure. Firstly, the hamiltonian L has a dou- 
bly degenerate ground eigenspace spanned by the vec- 
tors {|0)|Ojy), |l)|fljf)}, where |fiff) (respectively, \£Ik}) 
is the ground state of H (respectively, K). The hamil- 
tonian L then has a gap AE which is larger than some 
constant, irrespective of the number 77 of spins. 

To construct our final auxiliary system K, we consider 
the parameter-dependent hamiltonian 

M{6) = L + K V C >(0)<8>I C (4) 

where 

V{6) = (sin(0)|O)-cos(0)|l))(sin(6')(O|-cos(e')(l|). (5) 

We note that V(0) is positive semidefinite and it acts 
nontrivially only on the auxiliary spin. Thus, to find the 
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ground state of M(9) we can restrict our attention to 
subspace spanned by {|0)|Qjj), |1)|Qk-)}. 

The addition of the operator V(9) will perturb all of 
the eigenvectors of L. However, by Weyl's perturbation 
Theorem as long as K < the operator V(9) will 
not mix the ground subspace of L with the remaining 
eigenvectors of M(9); this subspace will always be sepa- 
rated by a gap from the rest of the spectrum. Further- 
more, the subspace itself is unchanged: only the eigen- 
vectors within this subspace change under the addition 
of V(6). Thus we fix k = ^f-, so M{9) = L + ^f-V{9). 
The matrix elements of V{9) in the ground eigenspace 
of L are given by V jk (6) = (^|F(Wfc>, j,k € {0,1}, 
where = \0)\Q H ) and |V>i> = (1)1^): 



sin 2 (9) 
— sin(#) cos(8)x 



sm(9) cos(9)x 
cos 2 (9) 



(6) 



The two eigenvalues and eigenvectors of this matrix cor- 
respond to the ground state and first excited state of 
M(9). The corresponding gap of V(9) is 



A(0) = ^/l-4sin 2 (0)cos 2 (0)(l- |af), (7) 

which has a minimum value equal to |x| at 9 — j, where 
x is defined by Eq. |T|). 

We think of the V(9) contribution in M{9) as polaris- 
ing the system L so it has a unique gapped ground state 



\n M (0)) = cos(0)|O}|fi H }+siii(0)|l}|fljO 



(8) 



By the discussion in the previous paragraph the gap 
above this ground state is always larger than AE J X \ . 

The idea behind our proof is now simple to state. We 
begin with the system K. in the ground state |Om(0)) = 
|0)c" \Qa) Mb) of Af(0) and adiabatically vary 9 from 
to j. The resulting ground state is |l)c|0#). This is a 
product state between the original chain C and C". We 
then discard the ancilla spin C to obtain the ground state 
of C. We approximate this exact adiabatic evolution with 
a unitary operator which acts nontrivially only on a small 
set A of spins across the boundary between A and B. 



IV. APPROXIMATING THE ADIABATIC 
DYNAMICS 

The adiabatic evolution of the ground state \Qm{9)) = 
W(0;O)|O M (O)) of M{9) is generated by the solution of 
the differential equation 



d9 



U(9;0) = [fi' M (0),n M (0)]W(0;O), 



(9) 



where £Im(S) = \^m(9)} (^m{9)\- This is an example of 
an exact adiabatic evolution (see App. [X] for more dis- 
cussion of exact adiabatic evolution). 

We approximate the exact adiabatic evolution U (9) by 
quasi- adiabatic evolution which, for us, is defined by the 



solution of the differential equation 



ds 



V(s;0) = 
X 7 (*) 



r M( s ; 



f dM(t 
\ ds 



where 



e 2 -< 



27T7 ' 



du )dtV{s;0), (10) 



(11) 



and 7 is some constant to be set later, and r t M (AT) = 
e ztM^Y e -jtAf ^g ee App. [5] for a discussion of quasi- 
adiabatic evolution and related adiabatic-like evolu- 
tions.) 

The infinitesimal generator 



k(s) = 



x 7 W 



u \ ds 



du dt (12) 



of quasi-adiabatic evolution is an operator which is ap- 
proximately local in a region Aj = {I \ d(l, C) < j} of 
2j + 1 sites surrounding the boundary between A and 
B. The way to see this intuitively is to recognise that 
d Q s s ^ is strictly local in the region Ai, and the operator 
Tu' 1 ^ \ ~§^ j i s a local operator evolved according to a 
local hamiltonian for some time which is approximately 
less than some constant 7. We make this intuition pre- 
cise by applying a Lieb-Robinson bound [3, [H, [TtJ 
(see [l8[ for a simple direct proof). The Lieb-Robinson 
bound reads (for a system with hamiltonian H) 



|[7f(A),iJ]||<|y| e - 



-v d(x,Y) / /h: 1 1 j 



(e^'-l), (13) 



for any two norm-1 operators A acting on site x and B 
acting on a subset Y of sites, with {x} n Y — 0, which 
are separated by a distance d(x, Y). The constants v and 
k are independent of n and depend only on \\h\\, which 
is an O(l) constant. 
What we do is define 



k a (s) 



f MA„( s ) f dM{s) 
\ ds 



(14) 



where 



^W ( . )= j" x ^ t )( l'r^ s \-)du)dt, (15) 



with 
and 



ith M A „(s) = I& ® H Aa + ^V(s) + ® Hi 

jeA a 

where hj = T J (h). Obviously k a (s) has support 
supp(fc a (s)) = A Q . 
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We want to show that the quasi-adiabatic dynamics 
generated by k a (s): 

d 



so that 7 > \ X \ \ E is sufficient to ensure that the approx- 
imation is exponentially small (in 7). 



ds 



V Aa (s;0)=ik a (s)V Aa (s;0), 



(17) 



With the choice a > 7 > 



TAE 



we find that the width 



[i] 



acts on 



are close to the quasi-adiabatic dynamics V(t;0) gener- 
ated by k(s). We do this by exploiting the inequality 

f \t\ 

||V(t;0)-V Ao (t;0)|| < / \\k(s) - k a (s)\\ds, 
Jo 

which is proved, for example, by exploiting the Lie- 
Trotter expansion. 

We now show how the Lieb- Robinson bound provides 
an estimate on the decay of ||fc a (s)||. Consider 



of the region A a that the unitary operator V A 
is given by , y± E , where c is some constant. 

To conclude our discussion we now show how to itera- 
tively use the procedure we've described in the previous 
paragraphs to construct an approximation to the ground 
state |f2ij„) of H n which can be stored (as a FCS) with 
resources scaling as a polynomial in n. 

Our first step is to start with a system A\B\ which 
consists of two copies A\ and A 2 of the chain C m on to 
sites with to > 



k(s) - k a (s)\ 



is given by 



\x\AE- 



The hamiltonian for this system 



.M(s) 



< 2 



< 2 



< 



l*r(*)l 
lx 7 (*)l 



(to') - t 
fW(m') 



Ma„ (s) . 



(to')) du ) dt 



K = H 

J *■ m. - L - L 77 



Isi + Ia! ®H m - S m I. 



(21) 



M a (s) 



K) 



du I dt 



< 



< 



lx 7 (*)l 



min{2|jm'||,ca 2 A£;e K|tl| ^ Q }rf 
'■ ' r -'dl A I / \ Xl {t)\\t\dt 



dt 



"dt + AE 



27T7 



\t\dt 



This system has a ground state equal to |O m )|f2 m ). We 
then follow the procedure described above, namely, ad- 
joining an ancilla spin C" between the two blocks, and 
then applying the approximate adiabatic evolution V^j 

to yield an approximation vj^ |fi TO )|fi m ) to the ground 
state |fi 2m ). 

We next iterate this procedure: we use two copies 
of the approximation V^' |O m )|f2 TO ) to approximate the 
state |il 2m )|il 2m ). This is the ground state of a new sys- 
tem A2B2 whose hamiltonian is given by 



C 2 Q 2 
27 2 . 



K 2m = H 2m <S> Is, + Ia 3 ® #2m - 5 2m I. 



(22) 



where to' = — = g ( s ' , and in the first line we ap- 
plied the triangle inequality, in the second line we applied 
the Lieb- Robinson bound, and in the third line we've bro- 
ken the integral into two pieces and applied the different 
regimes of the Lieb-Robinson bound separately with c 
some constant and we've used the fact that ||to'|| < AE. 
Thus, by choosing c < v/k we see that \\k(s) — fc Q (s)|| is 
decaying exponentially fast in a for a > 7. 

Thus we learn that the quasi-adiabatic dynamics 
V(^;0) are exponentially close (in operator norm) to 

a unitary operator vj^ (f;0) which acts nontrivially on 
only the region A Q . 

In order to complete our discussion and show that the 
state Vj^(§;0)|0(0)) is close to \Q H ) = W(f ; 0)|n(0)) we 
apply the triangle inequality to bound 



^ oy the discussion above, we find that this hamiltonian 
is adiabatically connected to Hn m . Because we are using 

\Sl m )\VL m ) for If2 2ro ) we 



an approximation |f2 2j 



must account for the error that arises when we take two 
copies of our approximation: 



'2m/ 



<2| 
= 2| 

<2'- 



2m/ 



\9.' 7 



j m ))|^2m) 



2m/ v \ 

2j 2 \x\ 2 AE 2 

\x\AE 



2-yAEe- 
(23) 



in 



H 



-VH (|;0) |fi(0)>|| < \\\n H )-V (|;0) |fi(0))||- 

£(f;o)|n(o))|| (19) 



Finally, using this upper bound we can compute the 
error between |fi4 m ) and our approximation (obtained 
from approximating the adiabatic continuation from 
|n 2m )|n 2 m) to |0 4m )): 



^oj |fi(o)>- 

We use Eq. (|A29j) from the appendix to bound the first 
term and Eq. (fl8|l to bound the second term: 

,-27 2 |a;| 2 A£; 2 



|| \^4m 

where we've defined 



Vf> 2 J)|fi 2m )||<2e(7)+e( 7 ), (24) 



- 2 7 2 |a;| 2 AB 2 
la; I AS 



■^AEe' 



(25) 



n 



H 



■v£(f;o)|n(o))|| 



< 



\x\AE 



- jAEe~ 



(20) 



To obtain an approximation to the ground state |fi n ) of 
H n we iterate the above procedure |~log 2 (n/m)] times. 
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The error resulting from applying the above procedure is 
given by 

= " | i AF + 7A^e-^ . (26) 
to m \ \x\Ah / 

If we fix some prespecified error e and demand that our 
approximation satisfies |||On) — l^n)ll < 6 then we need 
that 

^ max {N^< c ' l0g (7)}< (27) 

where c and c' are constants that only depend on \\h\\. 
Hence, we learn that, similarly, the unitary operators 

(§) that we apply at each stage k act on a collec- 
tion of a > c"7 spins, with c" some constant which only 
depends on \\h\\. 

After applying the iterative procedure described above 
we end up with the following representation for our ap- 
proximation to |O n ): 

\Q n ) = w\n m )\n m )---\n m ) (28) 

where 

W = V Acv (1)Va„(2) ' ' ' VA„([log 2 (n/m)l-l), (29) 

with A a (J) = {k | d(k,mj) < a}. 

The representation Eq. (|28p is equivalent [lj| [2(| to a 
finitely correlated state requiring a number of degrees of 
freedom which scale as n2 c l A =l, where c is some constant. 
Alternatively, it is clear that the representation Eq. (j2"5|) 
is already in a form useful for extracting local proper- 
ties: the expectation values of local operators such as 
correlators are easy to compute using the representation 
Eq. (|28[) . Finally, it is worth noting that our represen- 
tation is also exactly in the form of a simple instance of 
the multiscale entanglement renormalisation ansatz in- 
troduced in plj . 

V. CONCLUSIONS AND FUTURE 
DIRECTIONS 

In this paper we have shown how a class of noncrit- 
ical ID quantum spin systems are adiabatically con- 
nected to a ID quantum spin system of — noninteract- 
ing quantum spins with local dimension 2 m . As long as 

to 3> max | y x ^ E ,c' log (7) j, where c and c' are con- 
stants which only depend on ||/i||, then the ground state 
|fi n ) of H n can be approximated efficiently. This re- 
sult bears a superficial resemblance to a naive applica- 
tion of real-space renormalisation: in real-space renor- 
malisation one argues that blocks of spins should be 
effectively noninteracting. (I.e., after log(-^g) renormali- 
sation group transformations we should be very close to a 
trivial fixed point.) While this intuition is clear it seems 



to be extremely challenging to put this intuition on a 
rigourous footing. 

One question presents itself at this point: for noncriti- 
cal spin systems with a gap AE do the boundary effects 
persist only a distance into the interior of the system? 
Despite the plausibility of this statement we've been un- 
able to prove it; there are many counterexample systems 
which appear to ruin the most obvious approaches. 



APPENDIX A: ERRORS IN APPROXIMATIONS 
OF ADIABATIC DYNAMICS 

The dynamics generated by the adiabatic evolution of a 
quantum system is an important paradigm in the study of 
quantum mechanics. One question which is particularly 
pertinent for the simulation of quantum systems is: can 
adiabatic evolution be simulated by the time-dependent 
dynamics of some (possibly different) quantum system? 
In the case of adiabatic evolution the answer is yes: one 
method is to just slowly turn on the interactions on a 
timescale which is small compared to the gap of the sys- 
tem. While this is an entirely satisfactory solution in 
and of itself one limitation lies in the error of this ap- 
proximation; when adiabatic dynamics is simulated via 
this method the error decreases as an inverse polynomial 
in the timescale T of the slowly changing dynamics [12] • 
Thus we are inspired to search for methods which pro- 
vide better error scaling. One such method is provided by 
quasi- adiabatic evolution [23| (see also [13]). In this case 
we find that the error scales exponentially as a function 
of the timescale T set by the gap of the system. In this 
appendix we show this exponential scaling. Our discus- 
sion is set in the wider context of "adiabatic" evolutions 
whose differential generators can be arbitrary. 

The framework we describe in this section to discuss 
quasi-adiabatic evolutions is similar to that introduced 
by [H and 



1. Quasi-adiabatic evolutions 

In this subsection we introduce a general framework 
to describe "adiabatic" -like evolutions for quantum sys- 
tems. 

We consider adiabatic quantum evolution generated by 
a parameter-dependent hamiltonian H(s) as s is varied 
adiabatically from s = to s = 1. Thus we would like to 
understand the ground state |0(s)) of H(s). We do this 
by setting up a differential equation for |0(s)): 

±\n(s))=p>(s)\n(s)), (Ai) 

as 

where P'(s) = aj(|n(s))(fi(s)|) and we've set phases [H 
so that (f2'(s)|s2(s)) = 0. Because P'(s) is not antiher- 
mitian the dynamics generated by this equation are not 
unitary. 
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There are at least two ways to set up differential equa- 
tions for |fi(s)) which do generate unitary d yna mics. The 
first is via exact adiabatic evolution (see |25l . [2(| for a 
rigourous discussion of rather general results about ex- 
act adiabatic evolution): 



ds 



\n(s))=-[p(s),p'(s)]\n(s)). 



(A2) 



Because of the gap condition on H(s), the "hamiltonian" 
[P(s), P '(s)] for this dynamics is given by first-order sta- 
tionary perturbation theory: 



The following formula for Q(s) may be verified by writing 
in its eigenbasis and exploiting the L2 unitarity of 



JtH 



the fourier transform: 



Q(s) = £ Xj(Ej(s) ~ fi( S ))|^-( S ))<^( S )| (A5) 



How close is Q(s) to P(s)7 We measure distance in 
operator norm: 



[p(s),p'(s)] = \n(s))(n(s)\ 



dH(s) 



ds Q(s)l - H(s) 
8H{s) 



Q(s)I-H{s) ds 



-|fi(s))(fi(*)|, (A3) 



where J7(s) is the ground-state energy of H(s), and 
we define w 1 „, ^ via the Moore-Penrose inverse: 

0(s)K— H (s) 

We now define an infinitesimal generator for a quan- 
tum evolution which is meant to simulate the adiabatic 
evolution. We begin by specifying a function X7W which 
is an even real function whose fourier transform Xj is de- 
caying rapidly outside some region [—7,7], and which is 
normalised so that £7(0) = 1- We use this function to 
create an operator Q(s) which is meant to approximate 
P(s): 



\\Q(s) - P(s)\\ = 



X-y(z-n(s)). (A6) 



sup 

z<Esp(H(s))\fl(s) 



This formula shows us that, in the case where H(s) has 
a gap AE(s) > A, Q(s) is close P(s) in operator norm 
as long as |x 7 (2:)| decays rapidly for \z\ > A. 

In the case that X7W is an even real function whose 
fourier transform X7 is C°°, has compact support in 
[—7,7], and is normalised so that x 7 (0) = 1 with 7 < A 
to ensure that only the ground state appears on the RHS 
of (|A4[) we can recover exact adiabatic dynamics: we first 
use the Duhamcl formula 



ds 



e itH ^ = 1 I e 



(t-u)H(s) dH ( S \ ™H(s) dz . 

ds 



Q(s) 



Xl {t)e- ltn ^e UH{s) dt, (A4) to rewrite jUJ: 

-OO 



ds 



|fi(a)) 



.dft(s) 
ds 



t Xl (t)dt\Q(s)) + 1 Xl {i)e 



V«W (^^-) ) • " - . -AT) 



r 



Using the fact that X7W is an even function of t and 
cancelling phases we obtain 



ds 



Ms)) 



X<y(t) 



ds 



(A8) 

By integrating this expression for ^|r2(s)) in the energy 
eigenbasis of H (s) and using the assumed gap structure 
one can find that this expression is equivalent to the usual 
expression obtained from first-order perturbation theory: 



of exact adiabatic evolution. Thus we find 

|Jl(a)) =Te 1 ti K(s)ds \n(Q)) = U(t; 0)|fi(0)) (All) 
dti^j dt\Q(s)). where we define 

U(t;0) =Te t ti K( - s '> ds . (A12) 



±\{l( S )) 
ds 

We use the form 
imal generator 



dH(s) 



Ms)). 



(A9) 



n(s)l-H(s) ds 

to deduce a hermitian infinites- 



K(s) 



*r(*) 



H(s) fdH(s) 



ds 



du dt (A10) 



It is now easy to guess a form for the infinitesimal genera- 
tor L x (s) which is meant to approximate exact adiabatic 
dynamics: we just allow the function X7( s ) m (| A10|) to 
be arbitrary: 



L x As) 



(A13) 

(We drop the subscript on L x ^(s) from now on.) The 
corresponding dynamics are obtained by integration: 



V(t;0) =Te i ti L{s)ds . 



(A14) 
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We call the dynamics V(t; 0) a quasi- adiabatic evolution. 

With the definition (|A14[) of quasi-adiabatic evolution 
we now define the (time-dependent) state | 

|*(t)>=V(t;0)|fi(0)>. (A15) 

The idea is that |$(i)) should be very close to |0(t)) as 
long as \\Q{s) — P(s)\\ is small. We make this rigourous 
in the next section. 



2. Error in quasi-adiabatic evolution 

In this section we study the error S(t) — |||$(t)) — 
|f2(t))|| between the exact ground state of H(s) and the 
state |$(i)) generated by quasi-adiabatic evolution. 



To get started on an upper bound for 5(t) we define 

\6(t)) = MO)) ~ V(t; Q)t|0(i)) = |O(0)) - V(0; *)!"(*)>■ 

(A16) 

It is easy to see, using the unitary invariance of || • ||, that 
5(t) = \\\S(t))l 

We use the fundamental theorem of calculus to write 
\S(t)) as 



\S(t))=j^ s (V( S ;0y\n(s)))ds. (A17) 

Now we use the definitions of 4-\Q(s)) and ^V(s;0) to 
find 



W)) = [ ("*>(., 0)tL W |n W ) + V(.s; 



|0(s)) ds 



(A18) 



Thus, using the triangle inequality and unitary invariance 
of II ■ II, we find 



*(*) = lll*(*)>ll < 



iL(s)\ft(s)) 



dH(s) 



Q{s)l-H(s) ds 



ds. 
(A19) 



To make some headway on this expression we derive an 
expression for L(s)\fl(s)): 

iL( s )\n(s)) = 

1 1"! Xl{t) Of T " {s) i^r) du ) dm{s)Y (A20) 

By integration, we find 



iL{s)\Vt{s)) = 



H{s) - n(s)I Q(s)l - H(s) ) ds 



dH(s) 



(A21) 



where we define, as before, QCTfjrgT^ y i a the Moore- 
Penrose inverse: ^(s)i-g(s) = 0> ana - we define the 

operator \^(H(s) — fi(s)I) using the holomorphic func- 
tion calculus [27|: 



3=0 

(A22) 



r 



Substituting this expression for iL(s)\fl(s)) into (|A19[) 
we find 



S(t)< 



%{H{s) - fi(a)I) dH(s) 



H(s)-n(s)I ds 
We rewrite this as 



\n( s )) 



ds. (A23) 



where 



and 



H(s) - Q(s)l 



|r/(s)) ds, 

(A24) 

(A25) 



2"-l 



(Recall we are defining ^withTO 
inverse.) 

We now define = Vugb^^Ms)) 
and 77* = sup se r t i T)(s) to rewrite (|A24|) as 



(A26) 

via the Moore-Penrose 



< 


dH(s) 




ds 



IIK*))II<»7 
We also define 

/(*) = 



x-y(H( s ) - n(s)i) 



H(s) - Q(s)l 

x 7 (H(s)-n( s )i) 



ds. (A27) 



H(s) - fl(s)l 

X 7 (z - Q(s)) 



(A28) 



sup ni \ 

z£ Sp (H(s))\n(s) z-il(s) 
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and set /* = sup sg [ t] f( s ) to obtain our final estimate 

|||<5(*)> II <W*- (A29) 

The quantities ry» and /» can be separately upper 
bounded. In the case of spin systems which are adia- 
batically evolving we can generally bound 77* by 77, < cn, 
where c is a constant and n is the number of spins. For 
/*: in the case that our system has a gap A we use the 
cutoff function 

_ t 2 

X7 (t) = 6 -^-, (A30) 

which has fourier transform 

X»=e- 2 ^ 2 . (A31) 



We can now put together the gap structure of the spec- 
trum of H (s) and the fourier transform of x 7 (i) to obtain 
the bound 

e -2 7 2 A 2 

/. < A ■ (A32) 

This expression is exponentially decaying in 7, indeed, 
7 > -j- is sufficient to make this upper bound for /* 
arbitrarily small. Putting these two bounds for 77* and 
/* together we find that, for adiabatically evolving spin 
systems, 7 > ^y/log(n) is sufficient to reduce the error 
S(t) until it is arbitrarily small. 

It is easy to extend our discussion to systems H(s) 
which have a continuous gapless spectrum as long as 
there is some bound on the growth of the density of states 
near the ground state [13, El • 
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